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Abstract 

We consider preprocessing a set S of n points in the plane that are in convex position into a data 
structure supporting queries of the following form: given a point q and a directed line I in the plane, 
report the point of S that is farthest from (or, alternatively, nearest to) the point q subject to being to 
the left of line £. We present two data structures for this problem. The hrst data structure uses 
space and preprocessing time, and answers queries in 0(2^^® logn) time. The second data structure uses 
O(nlog^n) space and polynomial preprocessing time, and answers queries in O(logn) time. These are 
the first solutions to the problem with O(logn) query time and o(n^) space. 

In the process of developing the second data structure, we develop a new representation of nearest- 
point and farthest-point Voronoi diagrams of points in convex position. This representation supports 
insertion of new points in counterclockwise order using only O(logn) amortized pointer changes, subject 
to supporting 0(log n)-time point-location queries, even though every such update may make 0(n) 
combinatorial changes to the Voronoi diagram. This data structure is the hrst demonstration that 
deterministically and incrementally constructed Voronoi diagrams can be maintained in o{n) pointer 
changes per operation while keeping 0(logn)-time point-location queries. 


1 Introduction 

Line simplification is an important problem in the area of digital cartography |(lro91l IDen98l IMS92| . Given 
a polygonal chain P, the goal is to compute a simpler polygonal chain Q that provides a good approximation 
to P. Many variants of this problem arise depending on how one defines simpler and how one defines good 
approximation. Almost all of the known methods of approximation compute distances between P and Q. 
Therefore, preprocessing P in order to quickly answer distance queries is a common subproblem to most line 
simplification algorithms. 

Of particular relevance to our work is a line simplification algorithm proposed by Daescu et al. jDMSWO^ . 
Given a polygonal chain P = (pi,p 2 , • ■ • ,Pn), they show how to compute a subchain P' = 
with ii = 1 and im = n, such that each segment [pi^Pi^j^.^] of P' is a good approximation of the subchain of 
P from pi^ to Pij^i ■ The amount of error is determined by the point of the subchain that is farthest from 
the line segment [pi^Pi^j^^]. To compute this approximation efficiently, the key subproblem they solve is the 
following: 
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Problem 1 (Halfplane Farthest-Point Qneries) . Preprocess n points pi, p 2 , • ■ •, Pn in convex position 
in the plane into a data structure supporting the following query: given a point q and a directed line i in the 
plane, report the point pi that is farthest from q subject to being to the left of line 1. 

Daescu et al. |DMSW()'H] show that, with 0(n log n) preprocessing time and space, these queries can be 
answered in O(log^n) time. On the other hand, a naive approach achieves O(logn) query time by using 
0(n^) preprocessing time and 0(n^) space. The open question they posed is whether O(logn) query time 
can be obtained with a data structure using subcubic and preferably subquadratic space. 

In this paper, we solve this problem with two data structures. The first, relatively simple data structure 
uses 0(n^+®) preprocessing time and space, and answers queries in 0(2^/® log n) time. The second, more 
sophisticated data structure uses 0(n log^ n) space and polynomial preprocessing time, and answers queries 
in 0(log n) time. Both of our data structures apply equally well to halfplane farthest-point queries, described 
above, as well as the opposite problem of halfplane nearest-point queries. Together we refer to these queries 
as halfplane proximity queries. 

Dynamic Voronoi diagrams. An independent contribution of the second data structure is that it pro¬ 
vides a new efficient representation for maintaining the nearest-point or farthest-point Voronoi diagram of 
a dynamic set of points. So far, point location in dynamic planar Voronoi diagrams has proved difficult 
because the complexity of the changes to the Voronoi diagram or Delaunay triangulation for an insertion can 
be linear at any one step. The randomized incremental construction avoids this worst-case behavior through 
randomization. However, for the deterministic insertion of points, the linear worst-case behavior cannot 
be avoided, even if the points being incrementally added are in convex position, and are added in order 
(say, counterclockwise). For this specific case, we give a representation of a (nearest-point or farthest-point) 
Voronoi diagram that supports 0(logn)-time point location in the diagram while requiring only O(logn) 
amortized pointer changes in the structure for each update. So as not to oversell this result, we note that we 
do not have an efficient method of determining which pointers to change (it takes 0(n) time per change), so 
the significance of this representation is that it serves as a proof of the existence of an encoding of Voronoi 
diagrams that can be modified with few changes to the encodings while still supporting point location queries. 
However, we believe that our combinatorial observations about Voronoi diagrams will help lead to efficient 
dynamic Voronoi diagrams with fast queries. 

Currently, the best incremental data structure supporting nearest-neighbor queries (one interpretation 
of “dynamic Voronoi diagrams”) supports queries and insertions in 0(log^ n/loglogn). This result uses 
techniques for decomposable search problems described by Overmars j()ve83| : see |CT92| . Recently, Chan 
|Cha06j developed a randomized data structure supporting nearest-neighbor queries in 0(log^ n) time, in¬ 
sertions in O(log^n) expected amortized time, and deletions in 0(log®n) expected amortized time. 

2 A Simple Data Structure 

In this section, we prove the following theorem: 

Theorem 2. There is a data structure for halfplane proximity queries on a static set of n points in convex 
position that achieves 0(2^/"^ log n) query time using 0(n^+®) space and preprocessing. 

Our proof is based on starting from the naive 0(n^)-space data structure mentioned in the introduction, 
and then repeatedly apply a space-reducing transformation. We assume that either all queries are halfplane 
farthest-point queries or all queries are halfplane nearest-point queries; otherwise, we can simply build two 
data structures, one for each type of query. 

Both the starting data structure and the reduction use Voronoi diagrams as their basic primitive. More 
precisely, we use the farthest-site Voronoi diagram for the case of halfplane farthest-point queries, and the 
nearest-site Voronoi diagram for the case of halfplane nearest-point queries. When the points are in convex 
position and given in counterclockwise order, Aggarwal et al. |ACSS89| showed that either Voronoi diagram 
can be constructed in linear time. Answering point-location queries in either Voronoi diagram of points in 
convex position can be done in O(logn) time using 0(n) preprocessing and space | E(fS8b| . 
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Lemma 3. There is a static data structure for halfplane proximity queries on a static set of n points in 
convex position, called Okey, that achieves O(logn) query time using 0{n^) space and preprocessing. 

Proof. Let pi,p 2 ,. ■ ■ ,Pn denote the n points in convex position in counterclockwise order. The Okey data 
structure consists of one Voronoi diagram V(i,j) for every subsequence pi,pi+i, ... ,pj of points, where 
indices are treated modulo n. The space and preprocessing is thus 0{n^). 

To answer a halfplane proximity query for a point q and a directed line I, we first find the subsequence of 
points on the left of line i. In O(logn) time, we can find the two edges {pi,pi+i) and (j)j,pj^i) of the convex 
hull that are intersected by the query line |0’R,98[ Section 7.9.1]. Then, depending on the orientation of the 
line i (i.e., which edge is struck first), we can decide between the two possible intervals: Pi+i,Pi+ 2 y ■ ■ ■ iPj or 
Pj+i,Pj+ 2 ,Pi- Then we locate q in the appropriate Voronoi diagram, either V{i + l,j) or V{j + l,i), and 
return the site pk that generated the corresponding Voronoi region. The total query time is O(logn). □ 

Transform 4. Given any static data structure V for halfplane proximity queries on a static set of n points 
in convex position that achieves Q{n) query time using M(n) space and preprocessing, and for any parameter 
m < n, there is a static data structure for halfplane proximity queries on a static set of n points in convex 
position, called V-Dokey, that achieves 2Q{n) + O(logn) query time using \n/m~\ M{m) + 0{n^/m) space 
and preprocessing. 

Proof. Let pi,p 2 ,. ■ ■ ,Pn be the n points in convex position in counterclockwise order. We define the [n/m] 
breakpoints to be the points pk with k = 1 (mod m), i.e., the points Pim+i for i € {0,1,..., [n/m] — 1}. 
The data structure consists of two substructures: 

T> Substructure: We construct an instance of the data structure V on the half-open interval of points 
between every consecutive pair of breakpoints. More precisely, for each i € {0,1,..., \njrrT\ — 1}, we 
construct an instance of T) on the points Pim+iTPim+ 2 , ■ ■ ■ ,Pniin{n,{i+i)m}- These structures require 
|'n/m](M(m) -I- 0(1)) space and preprocessing. 

Vorouoi Substructure: For each breakpoint pk, we construct Voronoi diagrams on all intervals of points of 
length an exact power of two with one endpoint at pk. More precisely, for each i G {0,1,..., [n/m] —1}, 
and for each j G {0,1,..., [lognj}, we construct the Voronoi diagram on the points pim+i+r for r G 
{0,1,..., 2-^ — 1}, and we construct the Voronoi diagram on the points pim+i-r for r G {0,1,..., 2-’ — 1}, 
where indices are treated modulo n. The space and preprocessing requirements for these Voronoi 
diagrams are 



Overall, the space and preprocessing required for 2?-Dokey is |"n/m] M{m) + 0(r\f jrri) as claimed. 

It remains to show how we can use X>-Dokey to answer halfplane proximity queries in 2 Q(ri) O(logn) 
time. Suppose that we are given a point q and a directed line i. As described in the proof of Lemma 01 
in O(logn) time, we can find the interval pi,pi+i,... ,pj of points to the left of line i. If this interval 
contains no breakpoints, then it is contained in the interval of a I? substructure, so we can answer the 
query in Q{n) time by asking the P substructure the same halfplane proximity query. Otherwise, let pii 
and pji be the first and the last breakpoints in the interval, respectively. We ask the P substructure 
immediately preceding pi> (representing the interval pii-rn,Pi'-m+i, ■ ■ ■ if i' > 0, or else the interval 

P(\n/ 7 ri\-i)m+iTP(\n/rn\-i)Ta+ 2 i ■ ■ ■,Pn) and the P Substructure immediately succeeding pji (representing the 
interval pj',pj'+i,...,Pmin{n,j’ +m-i}) the same halfplane proximity query. These queries cover the ranges 
Pi,Pi+i,... ,Pi>-i and pj>,pj'+i,... ,pj. To cover the remaining range pii,pi’+i,... ,pji between the two 
breakpoints, we use the property that any interval can be covered (with overlap) by two intervals of length 
an exact power of two. Namely, let k = 2 dsL )J^ where the difference j' — i' accounts for wraparound 
modulo n. We query q in the Voronoi diagram on the interval pf ,pi>+i,... ,pi>+k and in the Voronoi 
diagram on the interval pj>,pji-i,... ,pji-k- Together, the four queries cover (with overlap) the desired 
interval pi,pi^i,... ,pj. Among the four results from the four queries, we return the best (either farthest or 
nearest) relative to point q. □ 
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By starting with the data structure Okey of Lemma|21 and repeatedly applying the Dokey transformation 
of Transformational we obtain the structure Okey-Dokey-Dokey-Dokey-..., or Okey-Dokey^, which leads 
to the following: 

Corollary 5. For every integer k > 1, Okey-Dokey^~^ is a data structure for halfplane proximity queries 
on a static set of n points in convex position that achieves 0(2^1ogn) query time using 
space and preprocessing. 

Proof. The proof is by induction on k. In the base case fc = 1, we can use the Okey data structure from 
Lemma 13 because (2fc + l)/(2fc — 1) = 3. For fc > 1, assume by induction that we have a data structure that 
achieves query time at most c (2^'“^ — 1) logn using space and preprocessing at most Assume 

that the constant c is at least twice as large as the constants implicit in the O notation in Transform 0] We 
apply the Dokey transformation from Transform^to this data structure, substituting m = 

Thus, n/m = and n^jm = The resulting query time is at most 2c(2^'“^ — 

l)logn + (c/2) logn < c{2^ — l)logn as desired. The resulting space and preprocessing time is at most 
(n/m+ + (cl2)n^lm = c + l)n + (c/2) ^ ^^{2k-3)/{2k-i) 

cn+ (c/2) for sufficiently large n, as desired. □ 

The space and preprocessing time of Okey-Dokey^“^ according to Corollary[3can be written as . 

For any given e > 0, we choose k = 1/2 + 1/e:. Then the space and preprocessing time are 0(n^+^) and the 
query time is 0(2^/® logn), proving Theorem|21 


3 Grappa Trees 

Our faster data structure for halfplane proximity queries requires the manipulation of binary trees with 
a fixed topology determined by a Voronoi diagram. To support efficient manipulation of such trees, we 
introduce a data structure called grappa trees. This data structure is a modification of Sleator and Tarjan’s 
link-cut trees |ST83| that supports some unusual additional operations. 

Definition 6. Grappa trees solve the following data-structural problem: maintain a forest of rooted binary 
trees with specified topology subject to 

T = Make-Tree(u): Create a new tree T with a single vertex v (not previously in another tree). 

T = Link(t>, w, d, , TTir): Given a vertex v in some tree Ty and the root w of a different tree Ty,, add an 
edge {v,w) to make w a child ofv, merging Ty and Ty, into a new tree T. The value d G {i,r} specifies 
whether w becomes a left or a right child of v; such a child should not have existed previously. The new 
edge (v,w) is assigned a left mark of and a right mark ofmy 

(Ti,T 2 ) = Cut(u,w): Delete the existing edge {v,w), causing the tree T containing it to split into two trees, 
Ti and T 2 . Here one endpoint of {v,w) becomes the root of the tree Ti that does not contain the root 
ofT. 

Mark-Right-Spine(r, m): Set the right mark of every edge on the right spine of tree T (i.e., the edge from 
the root of T to its right child, and recursively such edges in the right subtree ofT) to the new mark m, 
overwriting the previous right marks of these edges. 

(e, to/, to/) = Oracle-Search(r, Oe)' Search for the edge e in tree T. The data structure can find e only via 
oracle queries.' given two incident edges {u,v) and {v,w) in T, the oracle Oe(u,v,w,mt,mr,ra'^,m'y) 
determines in constant time which of the subtrees of T — v contains x.^ (Note that edges {u,v) and 
(v, w) are considered to exist inT — v, even though one of their endpoints has been removed.) The data 
structure provides the oracle with the left mark mg and the right mark my of {u,v), as well as the left 
mark m'f^ and the right mark to/ of {v,w), and at the end, it returns the left mark to/ and the right 
mark to/ of the found edge e. 

^Given the number of arguments, it is tempting to refer to the oracle as 0{A, B, D, G, I, L, S), but we will resist that 
temptation. 
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Theorem 7. There exists an 0{n)-space constant-in-degree pointer-machine data structure that maintains 
a forest of grappa trees and supports each operation in O(logn) worst-case time per operation, where n is 
the total size of the trees affected by the operation. 

Proof. Our grappa-tree data structure is based on the worst-case version of the link-cut tree data structure of 
Sleator and Tarjan Section 5]. This data structure maintains a forest of fixed-topology trees subject 

to Make-Tree, Link, Cut, and several other operations, each in O(logn) worst-case time per operation, and 
using 0{n) space. The data structure represents each tree in the forest by decomposing it into a set of 
disjoint paths connected by tree edges. Each path is in turn represented by a biased binary tree whose nodes 
represent the vertices of the path, ordered in the biased tree according to the depth along the path. Thus, 
higher vertices in the path correspond to nodes farther left in the biased tree. (In fact, the structure in 
|ST83| puts the vertices of the path at the leaf nodes of the biased tree and uses the internal nodes of the 
biased tree to represent the edges of path, but it is easy to modify to use this view.) Thus, there is a bijection 
between edges of a path and edges in the biased tree representing that path. The link-cut tree structure 
for a binary tree T can therefore be seen as a tree R, the representation tree, in which each node has at 
most four children (up to two for the biased tree representing a path and up to two for interconnecting the 
paths/biased trees), the nodes of R correspond bijectively to vertices of T, and the edges of R correspond 
bijectively to edges of T. Note that the endpoints of an edge in R may differ from the endpoints of the 
corresponding edge in T, although the parent endpoint of the edge in R always corresponds to one of the 
endpoints of the corresponding edge in T (which endpoint depends on whether the edge is a left edge or a 
right edge in R). For an edge {v,w) between parent v and child w in R, if we restrict to the nodes in the 
subtree of R rooted at v, then the restricted sets of nodes on either side of {v, w) in R correspond exactly to 
the restricted sets of vertices on either side of the corresponding edge in T. Furthermore, by suitable biasing 
as described in [ST83|, R has height O(logn). 

We augment the representation tree R to enable marking as follows. Because our tree T has bounded 
degree (an assumption not made in |ST88j L we can also explicitly store T (the parent, left child, and right 
child of each vertex) and crosslink between corresponding nodes/vertices and corresponding edges in the 
two structures. To each edge of T we add a left-mark field and a right-mark field. These fields contain 
the last explicitly stored marks for the edge, and for edges connecting two different paths {nonpath edges), 
they are accurate, while for edges on a path, the right-mark field may become out-of-date. In particular, 
Link(u, w, m^, mr) sets the left-mark and right-mark fields of the created edge {v,w) to the specified values 
mi and mr, respectively. To each internal node of a biased tree in R, we also add a right-mark field, which 
may be blank. When nonblank, this field represents bulk right markings that should be (but have not yet 
been) applied to the entire subtree of the biased tree rooted at this node. Thus, the actual right mark of an 
edge e on a path in T is implicitly the first nonblank right-mark field of a node along the path from the root 
of the biased tree representing the path to the corresponding edge of e in R, if there is such a field, or else 
the right-mark field of the edge e itself. 

We can maintain this augmentation as the representation tree R changes. Because the definition of the 
augmented values is relative to individual biased trees, we care only about modifications to biased trees 
themselves, not about the modifications to the edges between different biased trees that form the entire 
representation tree R. The link-cut data structure modifies biased trees according to rotations, splits, and 
concatenations. We can modify the implementation of all of these operations to propagate the mark fields, 
at the cost of an extra constant factor, in such a way that preserves the implicit marks of all edges in T. The 
idea is to push down node marks judiciously: whenever any operation visits a node v in R with a nonblank 
right-mark field, copy that value to the right-mark fields of the edges in T corresponding to the up to two 
children edges of v in the biased tree, copy the value to the right-mark fields of the up to two children nodes 
of V in the biased tree, and then blank out the field in the node v itself. Because operations on link-cut trees 
always start at the root of R and traverse along paths down from there, any nodes involved in the operation 
will have already cleared their mark fields before they actually get used, so the marks on the corresponding 
edges in T will be up-to-date. 

To implement Mark-Right-Spine(T, m), we visit all biased trees that represent paths with edges along the 
right spine of T. We start at the biased tree at the root of R (which contains the node corresponding to the 
root vertex of T), and we walk to the leftmost edge in the biased tree that corresponds to a left edge in T. 
(This edge can be found by a simple augmentation to store which subtrees in a biased tree contain edges 
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corresponding to left edges in T.) The edges that precede e in the biased tree correspond to the portion of 
the right spine of T on the corresponding path. Some of these edges preceding e are on or adjacent to the 
path we just walked, and we set the right-mark fields of the corresponding right-spine edges in T explicitly. 
The remainder of the edges preceding e are in subtrees hanging off the left of the path we just walked, and 
we set the right-mark fields of the root nodes of these subtrees. The total number of markings is at most 
proportional to the length of our walk. Then we continue to the child biased tree that precedes e, if any, 
which represents the continuation of the right spine of T. (Again this child biased tree can be found by a 
simple augmentation to store which subtrees in a biased tree have a descendant child biased tree.) Because 
R has height O(logn), the entire length of the walk and thus the total number of markings is O(logn). 

Given a query oracle Oe and a tree T, we can perform Oracle-Search by a tree walk in R starting at the 
root. Upon visiting a node v with two children, we apply the oracle to the two children edges of v, which 
correspond to two incident edges of T sharing the vertex corresponding to v. The result of this oracle call 
tells us whether e is the parent edge of v or else which of the two children subtrees contains e. Here we use 
the correspondence between the sides of an edge in R and the sides of the corresponding edge in T, assuming 
that we have already narrowed our search to the subtree of R rooted at v. On the other hand, if we visit 
a node v with just one child, we apply the oracle to the unique child edge and the parent edge of v, which 
again correspond to two incident edges of T sharing the vertex corresponding to v. This oracle call tells us 
whether e is the parent edge of t; or e is in the child subtree. Because R has height O(log'n), Oracle-Search 
queries run in O(logn) worst-case time. Furthermore, any edge we visit during this traversal from the root 
will have its left-mark field and right-mark field up-to-date, because the ancestor nodes above the edge will 
have already been visited and thus their fields will have already been propagated and blanked, so we can 
provide the oracle with the left and right marks of each query edge. Similarly, once we locate the edge e, we 
know its left and right labels. □ 

4 Rightification of a Tree: Flarbs 

The fixed-topology binary search tree maintained by our faster data structure for halfplane proximity queries 
changes in a particular way as we add sites to a Voronoi diagram. We delay the specific connection for 
now, and instead define the way in which the tree changes: a tree restructuring operation called a “flarb”. 
Then we bound the work required to implement a sequence of n flarbs by showing that the total number of 
pointers changes (i.e., the total number of parent/left-child and parent/right-child relationships that change) 
is 0(nlogn). Thus, for the remainder of this section, we use the term cost to refer to (a constant factor 
times) the number of pointer changes required to implement a tree-restructuring operation, not the actual 
running time of the implementation. This bound on cost will enable us to implement a sequence of n flarbs 
via 0(n log n) link and cut operations, for a total of 0(n log^ n) time. 

The flarb operation is parameterized by an “anchored subtree” which it transforms into a “rightmost 
path”. An anchored subtree S oi a, binary search tree T is a connected subgraph S' of T that includes the 
root of T. A right-leaning path in a binary search tree T is a path monotonically descending through the 
tree levels, always proceeding from a node to its right child. A rightmost path in T is a right-leaning path 
that starts at the root of T. 

The flarb operation^ of an anchored subtree S of a binary search tree T is a transformation of T defined 
as follows; refer to Figure ^ First, we create a new root node r with no right child and whose left child 
subtree is the previous instance of T ; call the resulting binary search tree T'. We extend the anchored subtree 

5 of T to an anchored subtree S' of T' by adding r to S. Now we re-arrange S' into a rightmost path on 
the same set of nodes, while maintaining the binary search tree order (in-order traversal) of all nodes. The 
resulting binary search tree T" is the result of flarbing S in T. 

Theorem 8. A sequence of n flarb operations, starting from an empty tree, can be implemented at a cost 
of Oifogn) amortized pointer changes per flarb. 

^ “Flarb” is a clever abbreviation of a long technical term whose meaning we cannot reveal for reasons we cannot comment 
on at the moment, perhaps simply due to lack of space or of the aforementioned purported meaning. Note that this notion of 
flarb is different from that of ICal05l . 
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Figure 1: An example of a flarb. The anchored subtree is highlighted. 


Proof. We use the potential method of amortized analysis, with a potential function inspired by the analysis 
of splay trees |ST85| . For any node a; in a tree T, let w(x) be the modified weight of the subtree rooted at a;, 
which is the number of nodes in the subtree plus the number of null pointers in the tree. In other words, 
we add dummy nodes as leaves in place of each null pointer in T, for the purpose of computing subtree size. 
Define i.p{x) = Ig • Clearly \‘f{x)\ < lg(2n — 1), because the smallest possible subtree contains no 

real nodes and one dummy node, and the largest possible subtree contains n — 1 real nodes and n dummy 
nodes. The potential of a tree T with n nodes is 3>(T) = with the sum taken over the (real) nodes 

X in T. Therefore, |‘I’(r)| = O(nlogn) for any tree T. 

For the purposes of the analysis, we use the following heavy-path decomposition of the tree. The heavy 
path from a node continues recursively to its child with the largest subtree, and the heavy-path decomposition 
is the natural decomposition of the tree into maximal heavy paths. Edges on heavy paths are called heavy 
edges, while all other edges (connecting two heavy paths) are called light edges. 

Outline. To analyze a flarb in a binary search tree T, we decompose the transformation into a sequence 
of several steps, and analyze each step separately. 

First, the addition of the new root node r can be performed by changing a constant number of pointers 
in the tree. Because (p{r) = lg(2n — 1), the amortized cost of this operation is trivially O(logn). Thus, in 
the remainder of the proof, we focus on the actual restructuring of the resulting anchored subtree S' into a 
rightmost path, a process we call rightification. 

At all times during rightification, the nodes constituting the original anchored subtree S' continue to 
form an anchored subtree of the current binary search tree, and for simplicity of notation we continue to 
denote the current such anchored subtree as S'. 

To implement rightification, we first execute several simplifying steps of two types, called “zig” and 
“zag”,^ in no particular order. Each such step has zero amortized cost. Any number of such operations 
might need to be performed and we stop when neither can be applied. At this point, the anchored subtree S' 
has a particular form and we perform a final operation, called a “stretch”, at the cost of O(logn) amortized 
pointer changes. This bound, together with the observation that the potential drop over any sequence of 
operations is 0(n log n), gives the theorem. We now describe the details of zig-zagging and stretching. 

The zig. A zig is executed whenever a light left edge is part of the anchored subtree S'; see FigureEl The 
zig operation simply involves a right rotation on the edge in question. The actual cost of a zig is 0(1), which 
we set to be 1 to ease the analysis. 

To analyze the change in potential, let A and B denote the two children subtrees of the lower endpoint 
of the edge, and let C denote the right child subtree of the upper endpoint of the edge. We use the same 
letters to denote the modified weight of the subtree. Because the edge is light, A -|- i? -|- 1 < O. Then the 

^Unlike most terminology in this paper, these terms are used for no particular reason. Cf. footnote|^ 
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Figure 2: A zig: The thick edge belongs to the anchored subtree S' and is light. 





potential change is 

A$ 


Therefore, the amortized cost of a zig is zero, as claimed. 

The zag. A zag is performed whenever there exists, within the anchored subtree S', a path that goes left 
one edge, right zero or more edges, and then left again one edge; see Figure 13 The zag operation performs a 
constant number of pointer changes to re-arrange the path in question into a right-leaning path. The actual 
cost of a zag is 0(1), which we again set to be 1 to ease the analysis. 

We now argue that a zag reduces the potential by at least 1. First, notice that the contribution to the 
potential of parent nodes of the trees Bi, B 2 , ■ ■ ■, Bk decreases after the execution of the zag because, in each 
case, the left subtree remains the same while the right subtree grows. We will argue that the contribution 
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Figure 4: The form of the anchored subtree S' before the final stretch. The thick light edges are light, and 
the thick black edges are heavy. 


of the remaining nodes decreases by at least 1. Indeed, 
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as claimed. 

The final stretch. After all possible zigs and zags have been exhausted, "we claim that the anchored 
subtree S' must have the form sho'wn in Figure 01 Indeed, any tree that has no light left edge and no 
right-leaning path delimited by t'wo left edges must have this form. In particular, because the rightmost 
path in this tree must be light, its length is at most lg(2n -|- 1). 

The final stretch operation, which completes the flarb, simply converts this tree into a rightmost path by 
effectively concatenating the subsidiary right-leaning paths, incorporating them into the main path. Only 
0(log n) actual pointer changes are required. The potential does not increase because left subtrees of every 
node shrink and right subtrees grow, if they change at all. Therefore, the amortized cost of the stretch is 
indeed O(logn). 

This concludes the proof of the theorem. □ 
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5 Transformations 


In this section we show how flarbs and grappa trees come together with a little work to give us the main 
result. In the next two transformations, we focus on the farthest-point case, but the proofs apply equally 
well to nearest-point. 

Transform 9. Given a grappa tree data structure supporting each operation in O(logn) worst-case time, 
and given a data structure to incrementally maintain a tree created by n flarbs with O(logn) amortized 
pointer changes per flarb, we can construct an 0{nlog^ n)-space data structure that supports 0{logn)-time 
farthest-point queries on any prefix of a sequence of points in convex position in counterclockwise order. 

Proof. We construct an incremental data structure that supports 0(logn)-time farthest-point queries on 
the current sequence of points, {pi,p 2 ,.. ■ ,Pn), and supports appending a new point Pn-i-i to the sequence 
provided that this change maintains the invariant that the vertices remain in convex position and in counter¬ 
clockwise order. Thus the insertion order equals the index order and equals the counterclockwise traversal 
order of a convex polygon. The data structure runs on a pointer machine in which each node has bounded 
in-degree. Thus we can apply the partial-persistence transform of |DSST89] and obtain the ability to sup¬ 
port farthest-point queries on any prefix of the inserted points in O(logn) time. The space usage becomes 
proportional to the number of pointer changes during the insertions. 

We consider the ordered tree T formed by the finite segments of the farthest-point Voronoi diagram, 
ignoring their precise geometry; see Figure [3 More precisely, the farthest-point Voronoi diagram 
Section 6.3] divides the plane into n cells by classifying each point q in the plane according to which of 
pi,P 2 ,... ,Pn is the farthest from q. The farthest-point Delaunay triangulation is the dual of the 

farthest-point Voronoi diagram, i.e., it triangulates the convex polygon with vertices pi,P 2 , ■ ■ ■ ,Pn by con¬ 
necting two vertices whenever the corresponding Voronoi cells share an edge. We consider the dual tree T of 
this farthest-point Delaunay triangulation of the convex polygon, i.e., the dual graph excluding the infinite 
region exterior to the convex polygon. Each edge in this tree corresponds to (a nongeometric representation 
of) a finite edge of the farthest-point Voronoi diagram, which is the bisector of two of the points pi and pj 
that are adjacent in the Delaunay triangulation. Each node in the tree represents a vertex in the farthest- 
point Voronoi diagram, or equivalently a triangle in the farthest-point Delaunay triangulation, and therefore 
has degree d < 3, where any degree deficit corresponds to 3 — d infinite rays in the farthest-point Voronoi 
diagram not represented in the tree T. 

We can view the tree T as a binary search tree as follows. First, we root the tree at the node corresponding 
to the unique triangle in the Delaunay triangulation bounded by the edge connecting the first inserted point 
Pi and the most recently inserted point We view the infinite ray emanating from the Voronoi vertex as 
the “parent edge” of this root node, defining the notion of left child versus right child of a node according 
to the counterclockwise order around the Voronoi vertex. (Note that this order is the opposite of the order 
defined by the triangulation, so in Figure^ (right), we draw T in mirror image so that its geometric notions 
of left and right match that of the Voronoi diagram.) Second, we assign keys to nodes consistent with the 
in-order traversal. For each tree node corresponding to a Delaunay triangle with vertices Pi,Pj,Pk, where 
i < j < k, we assign a key of j. In other words, we assign the median of the three vertex labels of the 
Delaunay triangle to be the key of the corresponding tree node. 

One way to view this key assignment is as follows. If we imagine adding an infinite rays in place of each 
absent child in the tree, and add an infinite ray in place of the absent parent of the root (the dashed lines 
in Figure El right), matching the counterclockwise order around the Voronoi vertex, then we decompose the 
plane into regions corresponding to Voronoi regions, each of which corresponds to a single point pi. All of 
the nodes bounding pfs region correspond to triangles incident to pi. We assign the key i to the unique such 
node in T that is closest to the root of T, or equivalently the least common ancestor of such nodes, which 
is the inflection point between two descendant paths that bound the region. Two exceptions are i = 1 and 
i = n: the vertices incident to pi are those on the left spine of T, and the vertices incident to are those 
on the right spine of T. 

In this view, we also define the left mark of an edge to be the label of the region to the left of the edge, 
and similarly for the right mark. Thus, the two marks of an edge define the two points pi and pj whose 
bisector line contains the Voronoi edge. If an edge is the left edge of its parent node, then the edge’s right 
mark is simply the key of that parent, because the right edge of the parent creates an inflection point at the 
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Figure 5: Adding vertex vs in counterclockwise order. Top: Before. Bottom: After. Left: Farthest-point 
Voronoi diagram and its dual, the Delaunay triangulation. Right: Delaunay triangulation and its dual, the 
tree T with attached infinite rays drawn as dashed lines, drawn in mirror image so that geometric left versus 
right matches the order in the Voronoi diagram. The root vertex of T and its parent edge are emboldened. 


parent. Similarly, if an edge is the right edge of its parent node, then the edge’s left mark is the key of that 
parent. Intuitively, in either case, if we walk up from the edge on its “underside”, then we immediately find 
a local maximum in the region. On the other hand, in either case, the other mark of the edge is the key 
of the parent node of the deepest ancestor edge that has the opposite orientation (left versus right): this 
bending point is the first inflection point we encounter as we walk up the tree on the “top side” of the edge. 
The tree T is not balanced, so we use a grappa tree to represent it and the left and right marks of edges. 

Next we consider the effect of inserting a new point Pn+i- As in the standard incremental algorithm 
for Delaunay construction |dBvKOS99l Section 9.3], we view the changes to the farthest-point Delaunay 
triangulation as first adding a triangle pi,pmPn+i and then flipping a sequence of edges to restore the 
farthest-point Delaunay property. The key property of the edge-flipping process is that all flipped edges 
end up incident to the newly inserted point Pn+i- Therefore these changes can be interpreted in the tree as 
adding a new root node, whose left child is the previous root, and then choosing a collection of nodes to move 
to the right path of the new root. This collection of nodes induces a connected subtree because the triangles 
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involved in the flips form a connected set. (In particular, the flipping algorithm considers the neighbors of a 
triangle for flipping only if the triangle was already involved in a flip.) Thus, the changes correspond exactly 
to a flarb, with the flexibility of the flarb operation encompassing the various possibilities of which edges 
get flipped to maintain the farthest-point Delaunay property. Another way to view the addition of Pn+i is 
directly in the Voronoi diagram. The point Pn+i will capture the region i?„+i for which Pn+i is the farthest 
neighbor. The region i?„+i is a convex polygon. Outside the Voronoi diagram is unchanged, so all 

edges of the new Voronoi diagram are either bisectors of the same two points as before, or are edges of Rn+i- 
In T after the flarb, Rn+i corresponds to the right spine. 

Each pointer change during a flarb operation can be implemented with one cut and one link operation. 
Therefore the grappa tree implements the 0{n log n) total pointer updates from flarb operations in 0{n log^ n) 
total pointer updates. It remains to update the marks on the edges. By the incremental Voronoi/Delaunay 
view above, the only edges for which these marks might change are the edges incident to the new region 
i.e., the edges on the right spine. We update the right marks on all of these edges by calling Mark-Right- 
Spine(r, n -|- 1). The left mark of each edge on the right spine is simply the key of the parent node of the 
edge. During the execution of the flarb, various right paths were cut and pasted together with cuts and links 
to form the final right spine. The edges on the final right spine that were originally part of a right path in T 
already had a left mark equal to the key of their parent node. Any other edges on the final right spine were 
just added via links, so their left marks can be set accordingly by specifying the right mi argument to Link. 
Thus, the total number of pointer updates remains 0{n\og^ n). This concludes the space bound of the data 
structure. 

To support farthest-point queries, it suffices to build an oracle for the grappa tree’s Oracle-Search. 
Specifically, given two incident edges (u,v) and {v,w), the oracle must determine which subtree of T — u 
has the answer to the farthest-point query. Using the two marks on the two edges, two of which must 
be identical, we can determine the three vertices pi, pj, and pk of the Delaunay triangle corresponding to 
vertex v in T. The vertex of the Voronoi diagram corresponding to v lies at the intersection of the three 
perpendicular bisectors between these three vertices of the Delaunay triangle. We draw three rays from this 
Voronoi vertex to each of the three corners of the Delaunay triangle. These three rays divide the plane 
into three sectors, and the Voronoi regions corresponding to the nodes in each subtree oi T — v lie entirely 
in one of these sectors, with exactly one subtree per sector. In constant time, we can decide which of the 
three sectors contains the query point q. The farthest-point Voronoi region containing the query point q is 
guaranteed to be incident to the corresponding subtree, and therefore we obtain a suitable answer for the 
oracle query. At the end, Oracle-Search will narrow the search to a specific edge of T, meaning that the 
query point q is in one of the two Voronoi regions incident to the corresponding Voronoi edge. In constant 
time, using the two labels on that edge of the tree, we can determine which side of the bisector contains q, 
and therefore which Voronoi region contains q, i.e., which point pi is farthest from q. 

This concludes the proof of the theorem. □ 

Transform 10. Given an 0(nlog^ n)-space data structure that supports 0{logn)-time farthest-point queries 
on any prefix of a sequence of n points ordered in convex position in counterclockwise order, we can construct 
an 0(n log^ n)-space data structure that supports 0(logn)-time farthest-point-left-of-line queries on n points 
in convex position. 

Proof. Let pi,p 2 , ■.. ,Pn denote the n points in counterclockwise order. 

First we observe that, using the given prefix structure, we can also build an 0(n log^ n)-space data 
structure that supports 0(logn)-time farthest-point queries on any suffix of a sequence of n points ordered 
in convex position in counterclockwise order. We simply reflect the points about a fixed axis, reverse the 
order of the points, and build the prefix structure, and then apply the same reflection transformation to 
query points before giving it to the structure. 

Next we observe that, in O(logn) time, we can find the interval pi,pi^i,... ,Pj of points that are to the 
left of the query line. This algorithm is described in the proof of Lemma |21 

We build a collection of prefix and suffix data structures, and answer a query, via a divide-and-conquer 
recursion. The top level of the recursion is special because the sequence pi,p 2 ,... ,Pn is cyclic. In this case 
we build a prefix structure and a suffix structure on this list of points. These structures can be used to solve 
any query interval that contains either pi or or both. Namely, if interval contains exactly one of pi or p„, 
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then the interval is a prefix or suffix of pi,p 2 , ■ ■ ■ ,Pn- Otherwise, the interval is the union of a prefix and a 
suffix, so we can query both structures and return the farther of the two answers. 

At the general level of recursion, we have an interval ... ,Pj of points, i < j, and a guarantee 

that any query interval reaching this level of recursion is strictly contained within this interval (excluding 
both Pi and pj). At the top level of recursion, i = I and j = n and we know that the interval contains 
neither pi nor as required. Let m = [(i + j)/2j be the point midway between i and j. We construct a 
suffix data structure on the left half of points, pi,pi^i,... ,Pm, and a prefix data structure on the right half 
of the points, pm+i,Pm+ 2 j ■ ■ ■ ,Pj- As above, these data structures can be used to solve any query interval 
that contains either pm or Pm+i or both (and satisfies the assumption of being strictly contained within the 
interval pi,pi^i,... ,pj). Then we recursively build data structures in the left half and in the right half for 
query intervals that do contain neither pm nor Pm+i- In a query, we only need to recurse in one of the halves; 
we can decide which half overlaps the query interval in constant time by comparing m with the indices of 
the endpoints of the query interval. In the base case, j = i or j = i + 1 and there are no query intervals 
because of the strict containment, so there is nothing to do. 

The recurrence for query time is T{n) = T{n/2) + 0(1) = O(logn). The recurrence for space of the 
prefix and suffix data structures is T{n) = 2T(n/2) + 0(n log^ n) = 0(n log^ n). □ 

Combining Theorems[7|and|Hlwith Transforms HHand llOl we obtain the following main result of our paper: 

Corollary 11. There is an 0{n\o^ n)-space data structure that supports 0{logn)-time halfplane proximity 
queries on n points in convex position. 

We also mention the implication in the area of dynamic Voronoi diagrams, which follows from combining 
Theorems |7| and |S1 with Transform m 

Corollary 12. There is an 0{n)-space data structure for maintaining a nearest-point or farthest-point 
Voronoi diagram of a sequence of points in convex position in counterclockwise order. The data structure 
supports inserting a new point at the end of the sequence, subject to preserving the invariants of convex 
position and counterclockwise order, in O(logn) amortized pointer changes per insertion; and supports point- 
location queries in O(logn) worst-case time. 


6 Open Problems and Conjectures 

Several intriguing open problems remain open. One obvious question is whether the O(nlog^n) space of 
our second data structure can be improved while keeping the optimal O(logn) query time. One specific 
conjecture in this direction is the following: 

Conjecture 13. A sequence of n flarb operations, starting from an empty tree, can be implemented at a 
cost o/0(l) amortized pointer changes per flarb. 

We have no reason to believe that our O(logn) amortized bound is tight. Reducing the bound to 0(1) 
amortized would shave off a O(logn) factor from our space and preprocessing time. More importantly, 
it would increase our understanding of dynamic Voronoi diagrams, reducing the O(logn) amortized update 
time in Corollarv ll2l to 0(1) amortized. The potential function we use is inherently logarithmic; a completely 
new idea is needed here for further progress. 

On the issue of improving our understanding of dynamic Voronoi diagrams, we pose the following problem: 

Open Problem 14. Is there a data structure for maintaining a Voronoi diagram of a set of points in convex 
position that allows point to be inserted in log*^^^^ n time while supporting O(logn) point location queries? 

Here we relax the condition that the points be inserted in counterclockwise order, but maintain the 
restriction that they be in convex position. Although our potential function does not give the result, it is 
possible that a slight variation of it does. 

Finally, it would be interesting to improve the construction time in our second data structure, in particular 
so that it completely subsumes the first data structure: 
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Open Problem 15. Can the pointer changes caused by a flarb be found and implemented in o{n) time, 
preferably n time? 

We have not been able to fully transform our combinatorial observations about the number of pointer 
changes into an efficient algorithm, because we lack efficient methods for finding which pointers change. 
Solving this question would improve our construction time by almost a linear factor, and would provide a 
reasonably efficient dynamic Voronoi data structure for inserting points in convex position in counterclockwise 
order. 
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